Coulomb blockade in superconducting quantum point contacts 
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Amplitude of the Coulomb blockade oscillations is calculated for a single-mode Josephson junction 
with arbitrary electron transparency D. It is shown that the Coulomb blockade is suppressed in 
ballistic junctions with D — ► 1. The suppression is described quantitatively as the Landau-Zener 
transition in imaginary time. 



PACS numbers: 73.23.Hk, 74.80.Fp, 73.40. Gk 

Coulomb blockade phenomena in mesoscopic conduc- 
tors have been actively studied during the past few years 
They arise from the interplay of discreteness, in 
units of electron charge e, of electric charge Q of a small 
conductor, and tunneling into the conductor. Coulomb 
blockade requires for its existence "localization" of the 
charge Q, the condition that implies that the trans- 
parency D of the tunnel barriers isolating the conductor 
is small, D <C 1. In ballistic junctions with D — ► 1 the 
charge can move freely in and out of the conductor and 
both the charge quantization and associated Coulomb 
blockade are suppressed. Until now, full quantitative un- 
derstanding of such a suppression has been worked out 
only for nearly-ballistic single-mode junctions between 
two normal conductors It has been shown that 

for conductors with the quasi-continuous energy spec- 
trum, amplitude of the Coulomb blockade oscillations 
vanishes as the junction reflection coefficient approaches 
zero: R = 1 — D — ► 0. The aim of this work was to 
study this problem for superconducting junctions, where 
the situation appears to be different. Coulomb block- 
ade oscillations arise in this case || from the formation 
of Bloch bands in the Josephson potential U (<p) periodic 
in the Josephson phase difference ip. Since the ballis- 
tic junctions also have periodic Josephson potential, one 
could expect that Coulomb blockade exists even in the 
ballistic regime. It is shown below that this expecta- 
tion is incorrect and, similarly to the normal case, the 
Coulomb blockade is completely suppressed when R — ► 0. 
The suppression can be described quantitatively as the 
Landau-Zener transition in imaginary time. 

Coulomb blockade in superconducting junctions can be 
conveniently discussed as the quantum dynamics of the 
Josephson phase difference ip. The standard Hamilto- 
nian for quantum dynamics of <p (see e.g., JjJjTj) consists 
of the coupling energy H c (ip) of the junction electrodes, 
which in the case of low-transparency junctions reduces 
to a simple Josephson potential U (<p) , and the charging 
energy (Q — q) 2 /2C, where C is the junction capacitance, 
q is the charge injected into the junction from external 
circuit, and Q is the charge transferred through the junc- 
tion. Coulomb blockade manifests itself as periodic os- 
cillations of the junction characteristics as a function of 



the charge q with the period 2e. These oscillations can 
take place either in time |6|,|8|], when the junction is biased 
with a dc current I, and q — It, or as thermodynamic 
oscillations ||, if one of the junction electrodes is an iso- 
lated island and the charge q is induced on the junction 
capacitance by external gate voltage V g coupled through 
a gate capacitance C g : q = C g V g . In both situations, 
oscillation amplitude is the same and can be found from 
the junction free energy F(q). 

For a single-mode junction with large transparency, 
studied in this work, the coupling energy H c (ip) can be 
represented similarly to the normal case B as a sum of 
the energies Hl,r of electrons with momenta ±kp mov- 
ing forward and backward through the junction, and a 
potential V responsible for scattering between these two 
directions of propagation. The energy of the forward- 
moving electrons in a superconductor can be written in 
the standard matrix form: 



Hi 



A(x) 



-ihvpd/dx A(x) 
A*(x) ivpd/dx 



A, x < 0, 
Ae iv , x>0, 



(1) 



where "F^ = (ip L p ipLl) is the creation operator for quasi- 
particles with momentum kp , and vp is the Fermi veloc- 
ity. Hpi is given by the same expression with Vp — * — vp. 
The pair potential A(x) can be written in the step-like 
form (Q) under the assumption that the characteristic 
junction length d is much smaller than the superconduc- 
tor coherence length Jivf/A. 

Below we limit ourselves to the case of adiabatic phase 
dynamics assuming that all energies, including character- 
istic charging energy Eq = (2e) 2 /2C and temperature 
T are much smaller than A. The adiabatic condition 
Eq , T -C A allows us to make several simplification in 
the junction Hamiltonian. First, since in this case quasi- 
partilces are not excited in the junction electrodes, the 
charge Q is carried only by Cooper pairs. This means 
that Q can be expressed directly in terms of the Joseph- 
son phase difference ip p0| : Q — —2eid/dtp. Even more 
importantly, the energy spectrum of electrons moving in 
the contact can be determined in this regime treating ip 
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as stationary. The Hamiltonian Hp + Hp> is reduced then 
to a sum of the quasiparticle energies Skip) of the occu- 
pied states, and the total junction Hamiltonian can be 
written as 



2C { i dtp 



(2) 



The spectrum of eigenenergies Sk{<p) is found by solv- 
ing the Bogolyubov-de Gennes (BdG) equations with the 
pair potential A(x) (0). It consists of the continuum of 
states at energies outside the gap, |e| > A, and two dis- 
crete states in the gap: 



Combined with eq. (||), this equation means that when 
tp increases from to 2-7T, +kp states in the continuum 
move up in energy, so that precisely one state is removed 
from the lower half of the continuum, e < — A, and is 
added to the upper half, e > A. Together with the shift 
of the subgap states this means that the whole spectrum 
of +kp states shifts by one state up in energy. The spec- 
trum of — kp states shifts by one state down. The change 
in energies in this process is infinitesimal for all states be- 
sides the two states (||) which move across the gap and 
change their energy by 2 A. 



e ± {tp) = tAcos^/2, V ± (x) = y/£/2 



1 

-i<p/2 



^±ik F x— £\x\ 



(3) 



where £ = (A/hvp) sin tp/2. In all these expressions tp £ 
[0, 27r], and they should be continued periodically in tp 
beyond this interval. In the regime of classical phase 
dynamics, only the subgap states (|^) contribute to both 
the dc Josephson current and the ac current at 

low voltages and one could expect that only they 
are relevant for quantum phase dynamics. In fact, as we 
will see below, the continuous part of the spectrum plays 
important role in determining the effective potential for 
the evolution of tp. 

The subgap states merge with the continuum when 
tp = mod(27r). Equation (g) shows that as tp varies from 
to 2ir the state with momentum kp moves across the 
energy gap from the lower half of the continuum, e < —A, 
to the upper half, e > A, while the — kp state moves 
in the opposite direction. The states in the continuum 
also shift up or down in a similar fashion. To see this 
we calculate the variation of the density of states p(e) 
with evolution of tp using the Friedel sum rule (see, e.g., 
0H): 



c) 2 



d P (e) = 

dtp 2n dtpde 



lndetS(e) , 



(4) 



where S(e) is the scattering matrix for scattering off the 
discontinuity of the pair potential A(x) ([[]). Straightfor- 
ward solution of the BdG equations shows that for +kp 
states 



S(e) 



1 



M(l-e^) 
(1 - a 2 )e^ 



(1- 



') 



(5) 



where a(e) = sign(e)(|e| — (e 2 — A 2 ) 1 / 2 )/A is the ampli- 
tude of Andreev reflection from a superconductor. Scat- 
tering matrix for momentum —kp is obtained from (|^) 
by replacing tp with —tp. 
From eq. (0) we get 



^- [ dtp-^- IndctSYe) = 
2Wo dtp K> 



M>A, 

kl = A. 



(6) 




cp/27i 



FIG. 1. Two branches of the Josephson potential in a bal- 
listic junction with transparency D — 1: one that corresponds 
to the equilibrium occupation of Andreev states at tp — 
(s = 1), and another with equilibrium at tp — 2n (s = —1). 

Such a motion of the energy spectrum determines the 
effective potential for the dynamics of tp in the Hamilto- 
nian (j^). At tp = 0, when there are no state in the gap, 
the equilibrium occupation of the eigenstates implies that 
at I < A all the states with e < —A are filled, while 
those with e > A are empty. Since the adiabatic vari- 
ation of tp does not induce transitions between different 
quasiparticle states, the shift of the energy spectrum with 
these occupation probabilities gives rise to the following 
aperiodic potential for tp (Fig. 1): 



C%) = X>fc(¥>)=A(2m+(-l) 
m = int(| ip \ /2tt) . 



m-\-l 



"p/2), (7) 



The rise of the potential (f?j) with tp means that the 
phase can increase beyond the points tp = mod(27r) only 
at the expense of creating quasiparticles in the junction 
electrodes. In the case of classical Josephson dynamics, 
this process generates real quasiparticles and creates dis- 
sipative component of the Josephson current [p^[ . The 
energy relaxation restores then the 27r-periodicity of all 
the junction characteristics. It should be noted that the 
potential (0) for quantum phase dynamics can not be ob- 
tained if one takes into account only the subgap states 
Jl6| . It is also worth mentioning that the mechanism of 
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the spectrum shift creating the potential (Q) is very sim- 
ilar to the mechanism of the chiral anomaly in the ID 
quantum electrodynamics - see, e.g., ]T^] . 

An important consequence of aperiodicity of the po- 
tential (0) is complete suppression of the Coulomb block- 
ade oscillations in ballistic junctions. Since the Coulomb 
blockade in superconducting junctions results from the 
formation of Bloch bands in a periodic Josephson poten- 
tial, the aperiodicity of the potential obviously suppresses 
the Coulomb blockade. However, the periodic nature of 
the potential and Coulomb blockade are restored by finite 
reflection in the junction. Indeed the aperiodicity of the 
potential (0) is the result of the transfer of one occupied 
+kF states from the energy range e < —A to s > A and 
one empty — kp state in the opposite direction as phase 
evolves from to 2ir. The backscattering term V in the 
Hamiltonian (||) couples ±kp states and prevents such a 
transfer. One can see this by looking at the two subgap 
states (||) which in absence of backscattering cross the 
gap range in the course of ip evolution: the occupied ip + 
state which moves from e < —A to e > A, and the empty 
ip~ state moving from e > A to e < —A. At ip = tt, the 
energies of these states coincide and V couples them effec- 
tively. If this coupling is sufficiently strong, the occupied 
state which starts off at ip = as ip + turns into ip~ after 
passing the point ip ~ it and moves back into the energy 
range e < — A. Similarly, the empty state starting at 
ip = as 'ip~ turns into ip + and goes back to e > A. This 
means that at ip = 2n all states with e < — A remain 
occupied while those with e > A remain empty, as at 
ip = 0. In this way the backscattering couples the branch 
of the Josephson potential with no quasiparticles at 
ip = to the one with no quasiparticles at ip = 2ir (the 
same as (^|) but shifted along the ip axis by 2tt - see Fig. 
1), thus creating the periodic low-energy branch of the 
potential. 

Quantitatively, the backscattering term in the Hamil- 
tonian (||) is V = J dxll(x)p(x), where U(x) is the po- 
tential profile along the junction and 

p(x) = ]T ^i R a 3 y L , R + {^ R a^ L e 2lk ^ + h.c.) 

L,R 

is the operator of electron density. Here and below cr's 
denote the Pauli matrices. To evaluate the backscatter- 
ing term V in the basis of the subgap states (|J) we use the 
fact that the characteristic range of the potential U(x) is 
defined by the junction length d which was assumed to 
be much smaller than the coherence length Hvf/A. We 
find then that the only nonvanishing matrix elements are 
those that couple the two branches of the potential: 

\V\$ + ) = ir Asimp/2, (8) 

where r = ~iU(2kp)/hvp is the reflection amplitude of 
the junction Q, and U (2kp) is the Fourier component of 
the potential U(x). At small r, the backscattering term 



is relevant only in the vicinity of ip — it, where it 
reduces to irA. Then, the junction Hamiltonian (||) for 
ip G [0, 2tt] takes the following form in the basis of two 
branches of the potential: 




0.0 1 1 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 

R(A/E C )" 2 

FIG. 2. The probability amplitude w ( |ll| ) for the Joseph- 
son phase difference ip to stay in the low-energy branch of the 
Josephson potential in junctions with small reflection coeffi- 
cient R <C 1. 

The width of the Bloch bands and associated with it 
amplitude of the Coulomb blockade oscillations depend 
on the probability amplitude w of staying on the low- 
energy periodic branch of the potential in the Hamil- 
tonian (|^). This amplitude is controlled by the usual 
Landau-Zener transition, the same as in the case of clas- 
sical phase dynamics |Q. The only difference with the 
classical case is that now the transition should take place 
in the course of ip motion under the potential barrier, i.e. 
in "imaginary time" . Indeed, in the quasiclassical ap- 
proximation, the stationary Schrodinger equation with 
the Hamiltonian (||) and energy e ~ — A describing the 
evolution of ip near the level-crossing point ip ~ tt is: 

2(E c /A) 1/2 diJj s /dx = -sxip s /2 + VRip- s , (10) 

where x = ip — tt, and s = ±1 is the potential branch 
index (Fig. 1). In eq. © we removed the phase & of 
the coupling terms in the Hamiltonian (^J) by the sim- 
ple unitary transformation ip s — > e lsB / 2 ip s . Equations 
( |l0| ) are the imaginary-time version of the equations de- 
scribing the regular Landau-Zener transitions, and their 
solution is provided by the parabolic cylinder functions. 
From the asymptotes of these functions |l9| we find that 
the probability amplitude w for the state s = 1 starting 
at x — > — oo to reach the state s — — 1 at x — > oo is: 

(ii) 



(0) 
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The amplitude w is plotted in Fig. 2. It tends to 1 at R ^> 
(Ec/A) 1 / 2 , while w ~ (27TA) 1 / 2 at i? < (^c/A) 1 / 2 . 
Since the amplitude of the Coulomb blockade oscillations 
is proportional to w, eq. ( pi] ) shows that similarly to the 
normal junctions, in the superconducting case these os- 
cillation vanish asi? 1 / 2 ati?^0. 

The low-energy periodic branch of the potential in the 
Hamiltonian (||) coincides with the classical stationary 
Josephson potential which for arbitrary junction trans- 
parency D is 0@ : U ( ( P) = _A [ 1 ~ Dsin 2 (tp/2)} 1 / 2 . 
For D larger than the small ratio Ec/ A, the charactristic 
magnitude of the potential U(ip) is larger than Eq, and 
one can find the first few eigenenergies e n for tp motion 
in this potential using the quasiclassical wavefunctions 
away from the potential minima at tp — , 2ix and match- 
ing them to the oscillator wavefunctions in the vicinity of 
these points. Taking into account that the wavefunctions 
should be periodic: ip{ip + 2tt) = tp(tp), we find that each 
oscillator eigenenergy acquires a small correction — S n : 
e n = hujp(n + 1/2) — S n , where 



5 = AbDw 



E c 



2n 2 AD 



1/4 



cos—, (12) 

e 



b 2n ( AD\" /2 
S n = (-l) n S —r — 
1 ' n\ \2E C J 

Here us p = (EcAD/2h) 1 ^ 2 is the frequency of small oscil- 
lations around the potential minima, and q' = q—eO/n is 
the induced charge shifted by the phase of the backscat- 
tering coupling. The numerical factors a and b in eq. 
( |l2| ) can be expressed in terms of elliptic integrals of the 
first and third kind Jll|, and are plotted as functions 
of the transparency D in Fig. 3. At D -C 1, a = 2\[2 
and 6 = 4, while at D -> 1: a = 8(V2 - 1) + Rhx^R, 
b = 8(V2 - 1). The width 5 n of the Bloch bands de- 
creases gradually with increasing D at intermediate D's 
because of the exponential factor in eq. ( |l2| ) reflecting 
the increase of the Josephson potential, and then rapidly 
goes to zero at D — > 1 together with the probability am- 
plitude w Q). 

Summing the corrections 5 n (12) over n, we can find 



the g-dependent part of the junction free energy at finite 
temperatures (T -C A) on the order of huj p : 



F(q) = -S (q)(l - e- 



-^)exp{-6 2 (|g 



1/2 



,-hij p /T 



} 



(13) 

The free energy ( |l3| ) determines the amplitude of the 
Coulomb blockade oscillations, for instance, oscillations 
of the voltage across the junction: V(q) — dF{q)/dq. It 
should be possible to observe the D-dependence of the 
amplitude of the Coulomb blockade oscillations ([l3|) ex- 
perimentally either in the controllable atomic point con- 
tacts |2Q] or in the semiconductor/superconductor het- 
erostructures pjj . 



In summary, we have studied the Coulomb blockade 
oscillations in single-mode Josephson junctions with ar- 
bitrary electron transparency D in the adiabatic limit 
Ec "C A. It was shown that the amplitude of these 
oscillations decreases steadily with increasing D at inter- 
mediate D's and then is rapidly suppressed (on the scale 
(Ec/A) 1 / 2 ) at-D ~ 1. The rapid suppression is described 
quantitatively by the amplitude ( pi] ) of the Landau-Zener 
transition between two branches of the Josephson poten- 
tial. 




FIG. 3. Exponent a and the preexponential factor b in 
the amplitude of the Coulomb blockade oscillations dli| ) in 
a single-mode Josephson junction as functions of the junction 
transparency D. 
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